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Abstract 

Correlators of magnetic and electric field energy density are investigated for SU(iV c ) 
gauge theory at high temperatures T. At separations z < 2/T the correlators are shown 
to be dominated by a power-law behavior even for finite gluon screening masses. This 
continuum behavior is well approximated on current 4 x 16 3 ~lattices in the perturbative 
limit and leads to a considerable overestimate of screening masses deduced from fitting 
the lattice correlators with conventional exponential forms. The use of extended sources 
and sinks to enhance the signal improves the situation for screening masses m S> T 
but leads to a largely uncontrolled error for masses less than T. In fact, we show that 
recent lattice QCD data of Grossmann et al., from which a magnetic screening mass 
tum ~ 2.9 T was deduced, may even be consistent with a vanishing actual magnetic 
screening mass. 
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1 Introduction and Conclusions 



In a recent lattice QCD study of SU(3) pure gauge theory [|l] the plaquette-plaquette cor- 
relation functions were measured at high temperatures. Below the deconfmement transition 
temperature (T c ~ 200 MeV) the observed exponential fall-off of those correlators was 
used to determine the mass spectrum of glueballs. Above T c the two lowest masses de- 
duced from the decay of the correlators were interpreted as twice the electric and magnetic 
screening mass in the gluon plasma. The interesting result reported in [|| that we focus 
on is the claim that just above the critical temperature, T rs 1.5 T c , the magnetic screen- 
ing mass, tum/T 2.9 ± 0.2, is apparently twice as large as the electric screening mass, 
tue/T Ri 1.4 ± 0.2. Our interest in these results stems from the sensitivity of certain trans- 
port properties and hence signatures of quark-gluon plasmas produced in ultrarelativistic 
nuclear collisions [§, § to the ratio of these screening masses. In Ref. Q it was shown that 
the color relaxation time in a plasma is 

t c = [ZasTfimE/mM)}- 1 , (1) 

where f(x) ~ log(x) as x — » oo. This time scale controls the lifetime of collective color 
excitations in the plasma such as the plasmon. Also next-to-leading-order corrections to 
the electric screening mass depend self-consistently on the ratio mg/mM f|. In lowest order 
perturbation theory, m E /T = g[(N c + Nf /2)/3] 1 / 2 while ttim/T = 0. The infrared mag- 
netic singularities in higher orders are however expected j7| to produce a non-perturbative 
m Ai/T ~ g 2 . In the weak-coupling limit (T ^> T c ), m E /m M ~ 1/g ^> 1, and there is a natu- 
ral ordering of the screening scales, T 3> gT ^> g 2 T. However, at T m 1.5 T c , m^/mi/ Ri 1/2 
|jlj], and the measured screening masses imply a reversal of that perturbative ordering. Since 
the temperatures accessible via nuclear collisions are of the order of 2T C , this inverse order- 
ing of screening scales could have important phenomenological implications. For example, 
it implies a higher color conductivity which in turn could decrease the produced transverse 
energy and the dilepton rate ||. 

The screening masses are also of general interest as another test of the viability of the 
perturbative picture of the quark-gluon plasma. On the one hand, the lattice entropy den- 
sity |J, color electric screening mass |10[, baryon susceptibility |]TTJ, and spatial dependence 



of higher mass meson and baryon correlators []12| are all consistent with the perturbative 
picture of the quark-gluon plasma, even near T c . On the other hand, the temperature depen- 
dence of the energy density and pressure deviate strongly from the ideal Stephan-Boltzmann 
behavior up to several times T c ]TB"|, [TjL |T5|, the decay of the scalar and pseudoscalar meson 



correlation functions cannot be reconciled with free propagation of the gg-pair |12[ , and 
gg-wavefunctions behave as if the quarks were strongly bound [IE|. 

The recent measurement [0] of rrtM ~ 2m,E poses a further problem for the perturbative 
picture of the gluon plasma. How can static color magnetic fields be screened on such a 
small distance scale? In a magnetic screening volume 47r/(3m 3 M ) there is less than 1/2 gluon 
on the average in the perturbative picture. Non-perturbative mechanisms, such as a very 
dense color magnetic monopole gas, would have to exist to account for that large a magnetic 
screening mass. However, there appears to be no consistent effective QCD action that could 
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generate such a monopole gas [T7] . 

In this work we aim to clarify how the lattice results of |J above T c may translate 
into continuum physics. To this end we assume that the lattice spacing is small enough 
to allow the plaquette correlators to be related to energy density correlators. The latter 
can be decomposed into multi-gluon exchange diagrams. We furthermore assume that the 
dominant contribution is given by the exchange of two dressed gluons. According to standard 
arguments, this would give rise to a decay of the correlator governed by twice the gluon 
screening mass, in agreement with the interpretation the authors of |[J have provided for 
their results. We first investigate the two-gluon exchange diagram in the interaction-free 
case to provide better insight on the screened case to be discussed subsequently. For distances 
z <C 1/7", i-e., on momentum scales where temperature effects are negligible, the correlators 
behave ~ z~ 5 , as expected from dimensional arguments. For z > 1/2T this behavior is 
changed to a decay ~ z~ 4 due to finite temperature effects. Then we study the effect of self- 
energy insertions. We show that the structure of the electric correlator is considerably more 
complex than that of the magnetic since the former mixes the longitudinal and transverse 
polarization effects. In both cases, for z < 2/T the correlators are found to be dominated by 
the power-law behavior encountered in the interaction-free case and not by the exponential 
decay expected for gluons with a finite screening mass. 

We then calculate the correlators on a lattice in the weak-coupling limit and show that 
the continuum limit is well approximated on a 4 x 16 3 -lattice. Consequently, the power-law 
behavior is still prominent Jl8| and leads to an overestimate of the masses deduced from 
conventional exponential-type fit functions. The power-law behavior is caused by the point 
source and sink of the correlator. However, in fuzzed links [0 were used which have 
the effect of "smearing" out the sources and sinks and partially remove the power-law en- 
hancement at small distances. We introduce a model for fuzzed links which incorporates 
the smearing and show that the relative residual error in determining the screening masses 
indeed becomes small if the latter are large compared to T. However, for a vanishing screen- 
ing mass this fuzzing model still leads to an apparent finite mass. In fact, we show that the 
current lattice data |l| may even be consistent with a vanishing magnetic screening mass ! 
The screening masses obtained from energy density correlators are therefore inconclusive and 



require confirmation from an independent measurement of e.g. the gluon propagator |20 . 

The remaining sections are organized as follows. In Section 2 we review in more detail 
the results of Ref. |lj and relate them to our investigations. In Section 3 the energy density 
correlators are decomposed in terms of multi-gluon exchange. Section 4 comprises a discus- 
sion of the interaction-free case. In Section 5 we discuss the correlators for dressed gluons. 
Section 6 is devoted to the transcription of the perturbative continuum results onto the lat- 
tice. In Section 7 we investigate a model for fuzzed links and their effect on the extraction 
of screening masses. 

We use natural units h — c — 1 and, except for Section 6 and 7, the Minkowski metric 
giiv — diag(H-, — , — , — ). This implies = (— ir, x) and = (— zcj n ,k), with the gluonic 
Matsubara frequencies u Tl = 2imT. 
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2 Plaquette correlators in lattice QCD 



In the perturbative picture of a quark-gluon plasma as a weakly interacting system, qq— 
(or meson) correlators are essentially given by a single gg-loop, since the mutual interaction 
between the quark and the antiquark is screened by the hot environment. Also, a quark 



acquires a dynamical mass m q ~ gT [21]. The correlator is roughly proportional to the 
square of the (dressed) quark propagator P2]| . For (large) spatial distances, the decay of the 
latter is governed by the lowest fermionic "energy" [(vrT) 2 + m^] 1 ' 2 [tiT being the smallest 
fermionic Matsubara frequency). Since 7r 2 3> g 2 we may neglect the quark mass as compared 
to the Matsubara frequency |23| and the spatial meson correlator decays like exp{— 2tczT} 



24fl . After correcting for lattice artefacts, the lattice data for vector and axial vector meson 
are found to be in good agreement with this picture The scalar /pseudoscalar channel is, 
however, an important exception in that it only decays ~ exp{— nzT}. This points towards 
a strong residual interaction in this channel, see e.g. []25], p6| . 

However, as the study of gg-wavefunctions has revealed |16| , even in the vector channel 
the quark-ant iquark pair may be strongly bound. The reason why the corresponding binding 
energy does not show up in the mass extracted from the correlator is the same why we could 
neglect the quark mass in our above consideration: it is negligible as compared to 2nT 
which originates from the lowest fermionic Matsubara frequency |16[ . It is therefore difficult 



to extract either the residual interaction effects or dynamical quark masses from spatial 
gg-correlation functions |27| . 



In view of this, the study of glueballs is of special interest: since the lowest Matsubara 
frequency for gluons is zero, the measurement of glueball correlators is not impaired by 
the above mentioned effect. They may thus provide more detailed information on residual 
interactions and the gluon screening mass. In []1| the irreducible representations for the 
glueball states on the lattice were identified and their masses measured via the decay of the 



corresponding spatial correlation functions It was found that above T c the lowest mass 



state no longer corresponds to a glueball with definite quantum numbers. Rather, it belongs 
to the (connected) correlator between plaquette operators of the form 

Pt = Pox + P<dy , (2) 

where 



1 



Oi 



Tr 



UoUiU^uf 



2N C 

is the elementary temporal plaquette with two spatial links Ui and two temporal links Uq, 

Up = exp {igaAfj} , 

with Afj, = and the generators of SU(iV c ) G b , 6 = 1,..., N 2 — 1. The lattice spacing a 

is a dynamical quantity in the sense that it decreases with decreasing coupling g according 
to the QCD /^-function @. Since g decreases with increasing temperature, at sufficiently 



high T the lattice spacing should be small enough to permit an expansion of the plaquette 
operator in powers of a. To lowest (non-trivial) order [^9j in a, the plaquette operator @ 
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is therefore proportional to 

P t ~Ti[8l + £l] . (3) 

where £, = jF i0 = F b G b = E\ G b is the zth component of the electric field strength [T^ is the 
field-strength tensor). Thus, (fj) is proportional to the color electric field energy density. As 
will be shown in detail in Section 3, the field energy density correlator can be decomposed 
in terms of multi-gluon exchange diagrams, the simplest of which is a one-loop, two-gluon 
exchange diagram. In the perturbative picture of the gluon plasma, this diagram is expected 
to give the dominant contribution to the energy density correlator, since mutual interactions 
between gluons are screened by the hot environment. In Section 5 we will show that in this 
case the mass determining the exponential decay of the energy density correlator is twice 
the gluon screening mass. It was found in |I[ that the mass extracted from the correlator 
between the Pt is consistent with the mass extracted from the Polyakov-loop correlation 
function, which for T ~ 1.5 T c agrees well with the perturbative result 2mg. This suggests 
that indeed the perturbative picture of the gluon plasma is valid, i.e., that the one-loop, 
two-gluon exchange diagram constitutes the dominant contribution to the P t -correlator at 
T m 1.5 T c . In our following continuum analysis, this will be our main working assumption. 
However, we emphasize that there is no rigorous proof that the respective temperature 
T ps 1.5 T c is sufficiently high and, consequently, the lattice spacing a sufficiently small so 
that the energy density correlators accurately reflect the plaquette correlators, and that the 
dominant contribution to the energy density correlators is the one-loop, two-gluon exchange. 

The second lowest mass extracted from the spectrum of all glueball-like excitations above 
T c was identified to belong essentially to the connected correlator between elementary spatial 
plaquettes 

Ps = Pxy i (4) 

which, to lowest non-trivial order in a, are proportional to the magnetic field energy density 
in z-direction: 

Ps ~ Tr E\ , (5) 

B z = T xy = F b y G b = B b z G b . In analogy to the electric correlator, the authors of [l] inter- 
preted this screening mass as twice the non-perturbative magnetic screening mass tum- 



3 Energy density correlators in terms of multi— gluon 
exchange 

As indicated in the last section, if one assumes that the dominant contribution to the energy 
density correlators is given by two-gluon exchange, then this allows an interpretation of the 
mass determined from the decay of the correlators as twice the gluon screening mass. In 
order to clarify this, we decompose the energy density correlators into multi-gluon exchange 
diagrams. We start with the magnetic correlator which has the simpler structure. The 
magnetic point-to-point correlator between r = (— ir, x) and r' = (— ir', x') reads 

M(r,r') = (Tr B 2 z (r) Tr B 2 z (r')) 
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Ti[G a G b ] Ti[G c G d ] (FUr)F b (r)FC(r>)F*(r> 



(6) 

The subscript "c" at the ensemble average brackets indicates that only the connected part 
of the respective expression is to be taken into account. Each color trace yields a factor 1/2 
and a Kronecker delta in the color indices. Inserting the definition of the non-Abelian field 
strength 

F; u = d,A a u - d v Al - gf abc AlAl , 
we arrive at the following decomposition 



M(r, r') = M (0) (r, r') + g 2 M (1) (r, r') + g 4 M (2) (r, r') 



(7) 



with 



2 r 



d x A a (r) - d y A a x (r) ' d x ,A b Jr') - d y ,A b x {r') 



V 

acd faef 



' jaca j}0 



d x A b y (r 



d y A b x (r) 



+ 



Al(r')A^r')Al(r')A^r') 
^4(r')-^4(0] 2 ^W^(r)^(r)4(r; 



(9) 



M(°)(r, r') 
M (1) (r,r') 

I jacd jt> e f 

M m (r,r') A _ _ _ % ..................... , /c . 

In order to calculate the ensemble averages we Fourier-transform the gluon fields according 
to 



d x A a y (r) - d y Al{r)] A x (r)A*(r) \d x ,A b y (r>) - d y ,A b x (r')] A e x (r')Af(r'] 



1 



f acd f a 9 h f hef f Ui /^c ( r )^( r )^( r ) M( r ) (r ') A*{l>)Altf) K (O ) (10) 



AJ(r) 



^=E e 



'11: 



The normalization makes the Fourier amplitudes A^(k) dimensionless numbers. In the ther- 
modynamic limit we have as usual 

d 3 k 



E=EE-E^ 



fcO k 



(2n) 



(V 



00 



An explicit calculation of the n-gluon amplitudes in is possible if we assume that the 

average is taken over a Gaussian ensemble, i.e., with an action of the form 

S[A] = -^^W {V- l )Z{k,k') A b u (k>) . (12) 
/J k,k' 

Note that D contains a gauge fixing term. The factor 1/T 2 is conventional in order that 
V^l agrees with the standard definition of the gluon propagator. The assumption of a 
Gaussian ensemble leads to the factorisation of an arbitrary n-gluon amplitude into a sum 
over products of all possible pairings of two-gluon amplitudes, i.e., gluon propagators, 
via 



= E ( A h A h) ■ ( A i n -l A in) 



pairs 



(K(k)Al(k')) 



t 2 v; b u (k,k') . 



(13) 
(14) 
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In the first line we used a collective index for momentum, Minkowski-, and color-indices. 
We furthermore assume translation invariance and global color neutrality so that 



(4) 

k,-k' 



With the abbreviation T/VJ2k = Ik we obtain 



M (0) (r,0) 



-i(k+l)-r 



k.l 



hi V yy (k) - 2 k x k y V yx (k) + k 2 V xx (k) 



i 2 x v yy {i)-2i x i y v yx (i) + i 2 y v xx {i) 



(15) 



and similar expressions for M^ 1 ' and . Each term in eq. flT5] ) has the graphical interpre- 
tation shown in Figs, la and Id. We represent the origin O and the point r by "vertices" 
(without the usual factor g !) and gluon propagators between them by full lines (for free 
gluons, i.e., V = V Q in (0), we use simple full lines as in Fig. la, for dressed gluons, we 
insert a self-energy "bubble" as in Fig. Id). A cross on the line corresponds to one extra mo- 
mentum factor or partial derivative. There are three topologically distinct two-loop graphs 
contributing to AfW and two three-loop graphs contributing to M^ 2 \ Examples of these 
are illustrated in Figs, lb and lc for T> = T> . However, we do not consider those further, in 
accordance with our working assumption that the one-loop, two-gluon exchange represents 
the dominant contribution to the energy density correlators. 

Other diagrams which will not be considered in the following arise from non-Gaussian 
fluctuations, which we have neglected in the factorisation flTB|), see for instance Fig. le. This 
generic ladder diagram could lead to bound states. In the perturbative picture of the gluon 
plasma, such a diagram would be suppressed since the ladder interactions are screened by 
the hot environment. Nevertheless, in view of the analysis of hadronic wavefunctions |16| 



one cannot a priori rule out that these graphs are not important near T c . Then, however, 
the exponential decay of the correlators would be determined by the mass of the exchanged 
bound gluon pair rather than by twice the screening mass of the two dressed gluons of 
Fig. Id. This points to an intrinsic ambiguity in the interpretation of the energy density 
correlators. Fig. If will be discussed in Section 7. 
The electric field energy density correlator, 



E(r, r') 



Tr 



Tr 



4V) + £,V)]) 



(16) 



can be treated in a completely analogous way. As in the magnetic case, our subsequent 
discussion will be restricted to the one-loop term 



#(r, 0) 



N 2 



1 



-i(k+l)-r 



k.l 



[{k Q ) 2 V i] {k) - k°k j V i0 (k) - WDP^k) + k i k j V 00 (k) 



h3=x,y 



x [(l°) 2 V lj {l) - l°PV M {l) - ri°V 0j (l) + /TO 00 (/)] • (17) 



We finally project the correlators onto a state with u n = k x = k y = 0. This is achieved 
by integrating them over the transverse r — x — y-plane at r, yielding the spatial magnetic 
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correlator 

rl/T rL/2 rL/2 

M(z) = dr dx dy M(r, 0) , (18) 

J J-L/2 J-L/2 

and similarly for the spatial electric correlator. The integration over the transverse plane 
eliminates three of the momentum sums in ([15]) and (|17|) via the identity 

l/T rL/2 rL/2 E 2 

dr dx dy e~ tk - r = e lkZz — 5 k o fi 5 k * t0 6 k y >0 . (19) 

J-L/2 J-L/2 1 

Explicit expressions will be provided in the next sections when the form of the gluon propa- 
gator is specified. 



4 The interaction— free case 

In this section we discuss the spatial correlators for the case of vanishing coupling g = 
(indicated by subscript at the respective quantities). In this case, the one-loop graph of 
Fig. la is the only term contributing to the correlators and is manifestly gauge invariant. 
Furthermore, T>^ u in eqs. fllB] , |T7|) reduces to the free gluon propagator which in the covariant 
gauge is given by 

gill* fcjijfc* 

For the integrand in (|lfj) we need the expression 



^(k) = »—-a — . (20) 



kl Vf{k) - 2 k x k y Vf{k) + k 2 y Vr(k) = ^ , (21) 

where k\ = k 2 x + ky. Note that the gauge parameter drops out, as expected. We now insert 
(|T5|) into ( |i8|) and employ (|T9D. Furthermore, in the thermodynamic limit 



1 _ r 00 dk 1 
T 2-^, = / ~wz 



L kl J-oo 2vr 



and 31 



Ufi ' ik z z 1 _ 1 -Ez 



e ^ - ^ = e -*>* Z>G (22) 

2vr E 2 + k 2 z 2E K ' 

where we have introduced the abbreviation E 2 = u 2 + k\. We thus arrive at 

Mf W = ^ri:r%^^e--. (23) 
2 „ Jo 27T AE 2 



We remark that the exponential decay factor in ( P2|) and in the integrand of ( [23]) is generic 
for spatial correlators. The smallest exponential factor (and thus the dominant term at large 
distances z) is that for the lowest Matsubara mode and zero transverse momentum. 



S 



For further purposes it is convenient to first consider arbitrary (even) powers of transverse 
momentum in the nominator of the integrand in (|23|) and to write E 2 = fi 2 + k 2 L . In our case, 
fi 2 = u 2 , but in general, when there is also a mass term in the propagator, /x 2 = u 2 + m 2 . 
Then, the integral reads 

i fW± e"-v^i = 1 1 Ci (2*H) , (24) 



8vr 7o /i 2 + k\ 8tt (2*)^' 

where the functions 



( v 2 - u 2 )i 

CM)= l ^ y ^ ^ 



are related to the exponential integral. We will only need [32 



C (u) = -Ei(-u)=e- u (^-J(u)) , (26) 

C x (u) = e' u {l + u)-u 2 C {u) = e- U (l + u 2 J{u)) , (27) 

C 2 (u) = e~ u (6 + 6u + 3u 2 + u 3 ) - 2u 2 e~ u (l + u) + u A C (u) 

= e~ u (6 + 6m + u 2 - u A J{u)) , (28) 



where 



Note that 



oo 

t 



J <"> s * e J^TW ■ (29) 



Cj(u) ~ , foru>0, (30) 
Ci(0) = 1,C 2 (0) = 6, andC (w)~lnw for u -> . (31) 



In eq. ( p3j) , ft = u n , and thus 



A/f (°) /V 2 — 1 1 00 



T 5 v ' 16vr (2zT) 4 



iV 2 - 1 1 



16tt (2zTY 



C 2 (0) + 2J2 C 2 (AixnzT) . (32) 



n=l 



With eq. (ESI), the sum over n can be further simplified using identities for the geometrical 
progression and its derivatives. For the following discussion, however, the form fl3~2| ) is more 
convenient. 

Let us first estimate the behavior of the correlator for different regions of zT. For very 
small distances zT -C 1, one expects that temperature effects do not play a role and that the 
behavior of the correlator is solely determined by dimensional arguments, as it would also 
be the case at T = 0. Since M(r, 0) ~ r~ 8 , we conclude that due to the integration over the 
t — x — y-plane in (|I8"D M(z) must be ~ z~ 5 |]T%|| . Indeed, taking that term in (|2"8]), which 
is dominant for small zT, and employing the geometrical progression one easily shows that 
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the sum over Matsubara frequencies n ^ in (^) diverges as 1/zT, and hence Mq^/T 5 has 
the expected (,2T)~ 5 -singularity. 

At larger zT the contributions of the non-zero Matsubara modes decay exponentially 



according to (|30|). Therefore, by virtue of (|3T|), the first term in curly brackets (corresponding 
to n — 0) is the dominant term for large zT. This changes the small distance (zT)~ 5 - 
behavior to a decay proportional to (zT)~ 4 . Thus, finite temperature enhances long distance 
correlations by one power of the distance as compared to the T = 0-case (in our case, from 
z~~ 5 to z~~ A ). This is a generic feature for all correlation functions at finite temperature and 
is essentially caused by the replacement of the time integration by T times a Matsubara sum 
in finite-temperature field theory 



These expectations are confirmed in Fig. 2a where we have plotted the full result (f3!|) and 
the contribution of the zero mode as a function of zT. For small zT, the non-zero Matsubara 
modes (~ (zT)~ 5 for zT — > 0) dominate over the n = 0-term (dotted line) which is only 
~ (zT)~ 4 , but already for zT « 1/2, they become negligibly small due to the exponential 
suppression, and the decay of the correlator is ~ (zT)~ A . 

Note that since ^(0) is constant, cf. eq. (0), the decay of the spatial correlator in the 
interaction-free case is polynomial, ~ (zT)~ 4 for zT > 1/2. However, as will be discussed in 
the next section, in the case of screening interactions the gluon will acquire a mass m. As 
mentioned above, this has the effect that // = u n — > (cu 2 + m 2 ) 1 / 2 , and the n = 0-term in 
(j32|) becomes C 2 (2zm). Consequently, it also decays exponentially on account of (^). This 
exponential decay is superimposed on the polynomial (zT)~ 4 -decay and of course dominates 
for large zT. At distances where zT as well as zm are small, however, the scales set by 
the temperature or the mass, respectively, cannot play a role and the correlator must again 
exhibit a ;? _5 -singularity as expected from dimensional arguments. 



Employing eq. (|17D , the analogue of flT8| ) for the electric case, and eqs. (\ZQ, |22|), the 
electric correlator reads 

<o) f s _ N 2 C - 1 _ r dk ± k ± k{ + 2u>l& 



V ; 2 ^Wo 2tt AE 2 y 1 



and can therefore be partly expressed in terms of the magnetic correlator, 

M (0) 

[zT) = ~^{zT) (34) 

oo 

J2 (^nzT) 2 { d(AiTnzT) + (ArmzT) 2 C (AnnzT)} . 



F (P) M (o) 

N?-l 

+ 



4tt (2zT) 4 



n=l 



Again, the sum can be analytically performed using identities for the geometrical progression. 

In Fig. 2b we show the electric correlator in comparison to the magnetic correlator. As 
expected from ([$(]), the additional terms in the electric correlator ( |34"D are negligible for 
zT > 1 and only influence the domain of 0.1 < zT < 1. For very small zT < 0.1, however, 
the electric correlator becomes equal to the magnetic (and thus has the same dominant 
singularity ~ (zT)" 5 ). The reason is that electric and magnetic fields are indistinguishable 
in Euclidean field theory at T = or at distances where the scale set by the temperature 
becomes irrelevant, respectively. 
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5 Correlators in the case of dressed propagators 



In this section we study the effect of the hot environment on the one-loop energy density 
correlators. To this end we assume that the two gluon propagators connecting O and r 
acquire a self-energy modification, see Fig. Id. The "dressed" propagator reads (in covariant 
gauge) 

with the two-space projectors 

k^k v 

V>» = A"" + — — , (36) 
qk" = gf* - - t>v» i (37) 

K 

where A M1/ = g^ v — u^u", A; M = /\^ v k u . is the (collective) 4-velocity of the matter under 
consideration. 

We now consider the combination of k l and T> % i appearing in eq. QT5| ) in the rest frame 
of matter, = (1, 0, 0, 0). We find the simple gauge invariant result 

kl vyy(k) - 2 k x k y vy\k) + k\ v*\k) = = k \ , (38) 

y -k 2 + 7r T uj z n + kj_ + k 2 z + ttt 



which should be compared with (|21f). The only modification thus far is the appearance of 
the "gluon mass" ^/t^t- In the general case, irj- depends non-trivially on k^. This may 
be seen considering for instance the gauge-invariant part of the one-loop gluon self-energy 
H [3H| . Only for k^ = (0,k), i.e., for the lowest Matsubara frequency, the dependence is 



trivial: 7ir(0, k) = for all k, indicating the absence of screening of static magnetic fields in 
first order perturbation theory. Static magnetic fields are screened in higher order in g, but 
this is not perturbatively calculable. On the other hand, lattice calculations may provide 
this information. As mentioned in Section 2, the results of [Q] were interpreted exactly in the 
way that they represent a non-perturbative measurement of the magnetic screening mass 
m M = (lim k ^ 7r T (0,k)) 1 / 2 . 

Given TTj<(k), in the general case one would have to perform the Fourier-transformation 
(p2|), with E replaced by 



E T = yjul + k 2 ± + 7T T , (39) 

and the subsequent integration over transverse momenta, cf. eq. (^), numerically. However, 
since there is up to now no way to obtain an analytical expression for 7r^, we may as well 
make the simplifying assumption that tc t = const., in order to proceed analytically. In this 
case, we anticipate (see also the discussion in the previous section) that the only modification 
of the final result in comparison to the interaction-free case will be an exponential damping 
of the correlator due to the gluon mass -Jt^t- 
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Indeed, for constant ttt, the result ( |2"2"D for the Fourier-transformation still applies (with 
E replaced by Ex) and we are left with the integral over transverse momenta, which is - 
after our groundwork in the previous section - readily expressed in terms of C2: 

M(o) at 2 — 1 1 f 00 ^1 

-^-w = -f^ (^7r 1 C2(2zV ^ ] + 2 S C2i n) \ ' (40) 

where = (a> 2 + ttt) 1 ^ 2 - On dimensional grounds ttt = CT 2 , where C is dimensionless. In 
the absence of any other scale C is independent of T, and thus the right-hand side of ( 40|) 
is still a function of zT only (as indicated in the argument on the left-hand side). 

For the following, we may assume C ~ 1. Then, for finite n we have 47r 2 n 2 T 2 « 40 n 2 T 2 ^> 
7r T ~ T 2 , and consequently fl^ ~ 27rnT. Thus, as required from dimensional arguments, we 
recover the (zT)~ 5 -singularity for zT — > 0, i.e., on scales where neither the temperature nor 
the gluon mass ^/iff (which is assumed to be of the order of T anyway) are relevant. In addi- 
tion, the exponential decay of the n 7^ 0-terms for large zT is similar to the interaction-free 
case. Consequently, these terms can be neglected for zT > 1/2. The most important modifi- 
cation as compared to fl32|) occurs for n = 0. In the interaction-free case, we had ^(0) = 6, 
cf. (j3i"D, leading to a purely polynomial decay ~ (zT)~ 4 of the correlator fl32|) for large zT. 
Now, however, due to the appearance of the gluon mass y 7 ^^, C 2 (2zy/7iT) ~ exp{— 2zy/Wr} 
according to (|30|) , i.e., an exponential decay is superimposed onto the polynomial decay and, 
of course, dominates at large zT. The corresponding slope is given by 2 v /7Pr, i.e., since 
m M = y/Wr, by twice the magnetic gluon screening mass, as argued in Section 2. However, 
we expect that in the region of small zT the (zT) _4 -behavior and the (;?T)~ 5 -singularity 
still influence the correlator to a large extent. 

For the electric correlator the combination of and T>^ v occurring in eq. (fl7|) reads with 

(H) 

(k°) 2 V ij (k) - k°k j V i0 (k) - k l k°V 0j (k) + k'kW^ik) 
5 ij - k' l ki/k 2 2 -k 2 k l k j 
-k 2 + 7i T ~ -k 2 + n L "k 2 " 

E 2 +k 2 \E 2 + k 2 kl + k 2 {E 2 +k 2 E 2 +k 2 Jf' [) 

where we have defined E\ = u 2 + k\ + 717, in analogy to ([39]). Again, this result is gauge 
invariant. 

It is noteworthy that not only ir^ enters the expression for the electric correlator, but also 
tit- The reason is that there is a transverse gluon A 1 involved in the definition of E(r, 0). 
One might therefore suspect that the exponential decay of the electric correlator will be 
eventually dominated by the smaller polarization function which could be ttt, at least at 
high temperatures where g 1 and ttt ~ g A T 2 <C ttl ~ 9 2 T 2 . Nevertheless, this is not the 
case: in eq. (f|l|) the respective terms are always accompanied by a factor of u 2 , thus they 



only contribute for non-zero Matsubara frequencies. As we have seen, these terms are in 
general strongly suppressed for large z. Furthermore, note that due to the structure of the 
projectors V^ u and Q^ v ', there is also a denominator k\ + k 2 . 
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The electric correlator is given by inserting ( |4"1"D into Again, further progress is 

possible under the assumption that both ttt and ttl are constant. As mentioned above, to 
one-loop order this is only true for = (0,k): 7Pr(0, k) = and 71^(0, k) 
Fourier-transformation can be performed using (^) and |35| 



m. 



The 



dx 



1 



a e 



-bz 



— he 



, a, b, z > . 



-oo 2tt (a 2 + x 2 ) (b 2 + x 2 ) 2ab (a 2 - b 2 ) 
We rearrange the final result in terms proportional to the various exponential decay factors, 




e -(k ± +E T )z + 



(7T T - 7T L ) 



-2E T z 



k 2 e~ 2k ^ z 



(42) 



Here we have also introduced fir 



71, 



(^) 4 #f 



(uj: 



tti) i . The integration over transverse momenta 
can be done analytically, but the result is extremely unwieldy. In fact, most of the terms are 
rapidly suppressed for increasing z, since they are accompanied by terms ~ u 2 exp{— zfl^ ,L }. 
As we will see, for a qualitative discussion of the correlator it is sufficient to retain only the 
last terms in the second and the third line of eq. (f42|) , 



dfcifci k 2 e~ 2k±z 



o 



E 2 L 



{2zf 

1 



C 2 (0) , 



(43) 
(44) 



For small z, the Matsubara sum over the terms (^4j) provides the z _5 -singularity expected 
from dimensional arguments. For zT > 1 and finite n, the term ([44]) becomes irrelevant, 
because it is exponentially suppressed ~ exp{ — 2zQ^} « exp{— AimzT} on account of flSOp . 
For n = 0, however, (a) it is the only term in ( f42|) which does not vanish and (b) it exhibits 
the weakest exponential decay ~ exp{— 2zy/¥i} of all terms in fl4~2"|). The slope of this decay 
is twice the electric screening mass ijie = \Pkl- The term (|43|) is only polynomially decaying 
(~ z~ A ) since £2(0) = 6, and will thus eventually dominate the long-distance behavior of the 
correlator (|4"2"D , provided n T ^ n L (which appears as a pref actor to ( [43]) in (ft"2])). Note that 
this term arises from the Fourier-transformation of the denominator /c 2 ^ + k 2 in (^TJ) which 
in turn originates from the specific structure of the projectors in eq. (BSI). For zT > 1 the 
approximate form of the electric correlator is 



g°) 



zT) 



N 2 



1 



16tt (2zT) z 



C 2 (2zv^) + 12 (7r T -7r L ) 2 ^ 



n=l 



M4 
n / 



(45) 



In Fig. 3a we show the electric and magnetic correlator, normalized to T 5 , as a function 
of zT for 7T£ = (1.4 T) 2 , it? = (2.9 T) 2 in comparison to the free case. This choice for the 
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polarization functions is motivated by the results of Ref. JTJ, where it was found that the 
electric correlator decays with the rate 2mg e m 2.8 T and the magnetic correlator 

with the rate 2tbm = 2\f^T ~ 5.8 T. 

One observes that the power-law behavior found in the interaction-free case dominates 
both correlators up to zT « 2. For larger zT the exponential decay prevails for the magnetic 
while for the electric correlator an anomalous power-law decay is observed due to the second 
term in (|45|) . This is shown in detail in Fig. 3b, where the various contributions to the 
electric correlator are plotted separately. The precise value of zT at which the power-law 
term dominates over the exponentially decaying first term in fl45f) depends on ttt and itl- 
This anomaly is an artefact from our simplifying assumption that the polarization functions 
do not vanish for large k and non-zero u n . In general we expect that for large zT the electric 
correlator decays exponentially on the scale 2m^. 



6 Lattice correlators in the weak— coupling limit 

Due to the periodicity of the lattice, a simulation can follow the decay of any correlator 
only up to z/L = 1/2, where L = aN a is the spatial extension of the lattice. This point 
corresponds, for any lattice with N a = 4N T , to zT = 2. What is surprising is that in the 
lattice simulations of QXJ] there is no hint of the power-law behavior which we have found to 
dominate the correlators for zT < 2 in our above continuum analysis. As we will show in 
the next section, this could be due to the use of fuzzed links |36 . 



In this section we quantify effects arising from the discretization and the finite lattice 
size without distortion from additional signal enhancement prescriptions. For this purpose 



it is sufficient to only consider the magnetic correlator. The transcription of eq. (|T|) onto a 



finite lattice proceeds along the lines given in [3~7|. After switching to the Euclidean metric, 
the main difference between the continuum and the lattice formulation is (a) that the gluon 
fields live in the middle of a link 

A » = Y,z lka{x+ae » l2)a A >) > ( 46 ) 

i\ a i\ T a k 

where e M is the lattice unit vector pointing in //-direction, and (b) that partial derivatives 
are replaced by forward finite differences, 

^ + , (47) 

The momenta in (HBT) fulfill 



h = , ji = 0, Ngr — 1 , i = x,y,z 



i 2n i° n at 1 

ko = — 7f , Jo = 0, ...,N T - 1 
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The free lattice action closely resembles the continuum action flT2D, 

1 



S [A] 



2MN T a A 



with 



rn 



sm 



kpCL 



(48) 



(49) 



where we have introduced a mass term. On the one hand, this allows to regularize any 
singularities due to constant field configurations. Usually, these singularities are avoided by 
omitting the ko = k x = k y = k z -term in momentum sums [57]. We will see that this is 
not necessary here: even for m = such a term is finite in the correlator (see below). On 
the other hand, allowing for a finite m from the outset provides a simple way to include a 
screening mass in the final expression. In the following, we will also make use of the relation 

(50) 



Given these facts, it is straightforward to derive the magnetic correlator 
Mf(r,0) = (Tr Jl(r) Tr ^(0) 







|sin 2 






2 _ 



D ,yy(k) - 2 sin 









0,c 


2 


-( 


k x a 




ky Qj 


sin 


_ 2 _ 




2 



1 1 v 

N*N T a 4 a 2 % 6 



ik-, 



(51) 



D 0>xy (k) + sin 2 



1^0,xx(k) 



The main difference as compared to (|i~5|) is that k l is replaced by 2 sin[/cja/2], as usual for 
bosons on the lattice. Summing over the r— ,x— , and y— sites at constant z in order to 
project onto zero momentum and employing ( f49|) we obtain 



M,S 0) (z)(aiV r ) 



iV 2 -l N, 



i\r2 S I jy X/ 



N 2 

o ko,k x ,k y 



a k z 



jk z z sin 2 [fc x a/2] +sin 2 [^a/2] 
H 2 +4E P sin 2 [fc p a/2] 



(52) 



The Fourier-transformation can be done using the identity 



N„-l 



E 



cos[2tt j z z/aN a ] 



cosh [2^(^/2 - z/a)} 
sinh[2u>] sinh[iV CT u;] 



(53) 



which can be proven by contour integration. Here u = hi[u + ywHl] and 

2 



UJ 



/ moV 



+ ^ sin 2 

p=r,x,y 



kpCL 



The reason why there is no singularity in (|52D in the case m = due to constant field 
configurations (ko = k x = k y = k z = 0) is that the term sin 2 [k x a/2] + sm 2 [k y a/2] regulates 
this singularity. Nevertheless, the limit m — > is subtle: putting m = /row t/ie outset, 
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the zero momentum term is finite, while performing the limit m — > after calculating the 
correlator, this term vanishes. Since the theory is not well-defined in the first case, we will 
always use the second prescription when considering the case m = in the following. 

In the case of a dressed gluon propagator one would in principle have to use the lattice 
analogues for the decomposition (|35|) and for the polarization functions ht,l @]- However, 
for an estimate of lattice effects and the general behavior of the correlator on the lattice it is 
certainly sufficient to simply assume a finite constant mass m = y/Wr in the above expression 
This assumption is supported by a comparison of ( |5TD and the continuum analogue 



with (BS 



The resulting correlator is shown in Fig. 4a for = and — (2.9 T) 2 in comparison 
to the corresponding continuum expression. One observes that the lattice ultraviolett cut-off 
q-1 regulates the z -5 -singularity at the origin, but that otherwise, on the 4x 16 3 -lattice used 
in [|TJ, the continuum limit is quite well reproduced. This, however, implies that the power- 
law behavior still influences the decay of the correlator. This has important consequences 
for the extraction of a screening mass. A fit with a simple exponential exp{— fi(z)z}, or, on a 
finite lattice with a hyperbolic cosine cosh[fi(z)(aN a /2 — z)\ |J (due to the lattice symmetry 
centered at z/L = 1/2), will necessarily lead to local masses fi(z) which are larger than lypnr 
(or 2y / ?TL in the electric case). Thus, screening masses extracted from these local masses via 
the simple relation \i{z) = 2rriE ) M overestimate the actual values. 

To quantify this we extract local masses from the correlators via 

M$\z) co^{ Z ){aN a /2-z)} 
M$\z-a) cosh\i{z){aN a /2 - (z - a))} ' 1 j 

which is the same procedure as used in Ref. jl]]. If this prescription worked correctly, the local 
mass should be two times the gluon mass, i.e., n(z)a = 2 ma = 2^Jttt/T 2 /N t = 1.45 = /i rea «a 
(iV T = 4 and tct = (2.9 T) 2 ). However, as one observes in the first row of Table 1 (see also 
Fig. 5a), even at z/a = 6 the local mass overestimates the actual value by 35% ! 

However, our results are not yet directly comparable to that of J3J. First of all, in that 
work sources and sinks of the correlators are plaquettes with a minimum edge length a, while 
( |52f) is essentially a point-to-point correlator. In addition, fuzzed links were used to enhance 
the signal at large z |[L9| . This technique effectively enlargens sources and sinks to the size 
of several plaquettes and also suppresses contributions of higher excited states at small z 
3"9]| , thus having an effect similar to wall sources commonly employed in the calculation of 



hadronic correlation functions. In the next section we investigate a model for fuzzing and 
apply it to our perturbative correlators. 

To conclude this section we remark that excited states are in our picture essentially 
states of non-zero transverse momentum and Matsubara frequency in the gluon loop in 
Figs. la,d. The summation over these excited states is nothing but the integration over 
transverse momenta and summation over Matsubara frequencies in ( p3|) or (42|). In this 
sense, the common notion that higher excited states enhance the point-to-point correlator 
at small distances is just reflected by our observation that the integration over transverse 
momenta and the summation over Matsubara frequencies leads to polynomial enhancement 
for small z. 
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Table 1 
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4.09 
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2 
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3.83 
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1.26 


0.96 


0.77 













2.83 
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2.30 


1.69 


1.31 


1.07 
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0.74 
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1.31 


1.05 
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0.81 


0.71 
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1.80 


1.10 
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0.81 


0.74 


0.66 
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7 Extended correlators 

Fuzzed links are obtained by replacing each link by a sum of products of two or more links 
in the vicinity of the old link |19| |. This roughly doubles the size of plaquettes at each new 



level of fuzzing and thus smears out the sources and sinks of the plaquette correlators. 

In our continuum analysis, we model fuzzing by introducing form factors, $(Ar) and 
$'(Ar'), for the smeared sink and source, respectively, and study the extended correlator 



M(z) = J 



J d 3 ri d 3 r 2 d 3 r[ d 3 r 2 $(Ar) $'(Ar') M{r 1 ,r 2 ,r[,r' 2 ) , (55) 
M(n,r 2) riy 2 ) = (Tt[B z (r l )B z (r 2 )]Tr[B z (r[)B z (r' 2 )]) c , (56) 



where Ar = r\ — r 2 , Ar' = r[ — r' 2 (Az = Az' = 0) and d 3 rj = dr^d^d?/,, i = 1,2, 
and similarly for d 3 r^. Note that, analogous to fuzzing [1], I9|| , smearing is limited to the 



transverse and temporal directions with z[ = z' 2 = and z\ = z 2 = z held fixed. Equation 
(|5|) reduces to ([18|) for the choice $(Ar) = 5 (3) (Ar), $'(Ar') = <5 (3) (Ar'). The prescription 
(^5|) is chosen such that (a) the only modification of the one-loop contribution (|23|) is the 
additonal appearance of the Fourier-transformed form factors $(— ui n , — kj_) and ^'(uv, kj_) 
and (b) its gauge invariance is preserved due to (|21~D. 

It is obvious that with suitable form factors one can suppress the contribution of large 
k± and n. Thus, our model correctly simulates the effect of fuzzing in that the small-z 
enhancement of the correlation function due to higher excited states can be reduced, which 
leaves the exponential decay less distorted. Diagramatically, the above prescription replaces 
Fig. Id with point source and sink with Fig. If, where the two dressed gluons are produced 
and absorbed from an extended source and sink, respectively. 

For example, wall sources correspond to the special choice $'(Ar') = T/L 2 = const., 
$(Ar) = 5^(Ar). This corresponds to $'(o; n ,kj_) = <5 n)0 > ^ = 1, which has the 
maximum effect of suppressing all higher transverse momenta and Matsubara frequencies. 
Unfortunately, wall sources are not applicable in our case, since the contribution of zero 
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transverse momentum has zero measure in the integral in ( |2"3"D and would thus lead to a 
trivial result for the correlator. 

We now investigate the influence of smeared sources and sinks on the lattice correlator 
(|52|). To this end, eq. (|55| ) is transcribed onto the lattice using standard discretization 
prescriptions. In |]J the sources and sinks prior to fuzzing are plaquettes of edge length a. 
Fuzzing links up to fuzzing level / increases the average edge length of the plaquette source 
and sink to 2 l a |19 |. Consequently, the sources and sinks have on the average the volume 



L T a x (L a a) 2 where L T = 2 1t , L a = 2 la . Thus, we choose as form factor 

Hr) = a3(Lr + 1)(Lff + 1)2 E o n E =q Knra 5 x , nxa 5 y , nya , (57) 

and the same for This corresponds to a homogeneous smearing of the source and sink over 
a cube of volume L T a x (L a a) 2 . Clearly many other possibilities for form factors are possible, 
but (p7|) is the most simple choice and will be adequate to prove our point, namely, that 
the extracted screening masses are extremely sensitive to the chosen smearing prescription 
and that smearing improves the extraction of screening masses only in the case when those 
are large compared to T. For small screening masses smearing will be seen to introduce 
uncontrolled errors. 

The fuzzing level used in is l T)(J = log 2 iV ri(T — 1, i.e., l a = 3 and l T = 1 for the 
4 x 16 3 -lattice under consideration. Consequently, L a = 8 and L T = 2. However, in order 
to study systematically the influence of our "fuzzing" prescription, we vary L a from to 
8 and L T from to 2 (see Table 1 and Figs. 4b, c, and 5a). Note that since sources and 
sinks of the correlators studied in |IJ are elementary plaquettes of edge length a prior to 
fuzzing, for instance an "unfuzzed" P s -correlator should actually not be compared to our 
magnetic point-to-point correlator of Section 6, but rather to an extended correlator with 
L T = 0, Lfj — 1. The Fourier-transform of (0) is 

$(cj„,k ± ) = a 3 e~ i{uJnT+k * x+k v y) $(r) 

r,x,y 

_ e -ia(uj n n T +k x n x +k y n y ) (58) 

(L T + 1)(L CT + l) 2 nr=0 nx>ny=0 



The result (|5*2] ) is modified in that the modulus squared of ( |5*%D appears additionally under 
the sum over ko, k x , k y . Using f40 |, one obtains the simple analytical expression 



,?/ i x, 2 _ ( sm[au n (L T + l)/2] sm[ak x (L a + l)/2] sm[ak y (L a + l)/2] \ 

1 [Un ' ±)l " \(L T + l)sm[au n /2] (L a + l)sm[ak x /2] (L c + 1) sm[ak y /2]J ' 1 j 

The correlator is displayed in Fig. 4b for ttt = (2.9 T) 2 and in Fig. 4c for ttt = for different 
degrees of smearing. One clearly observes a suppression of the polynomial enhancement at 
small z, which is proportional to the degree of smearing of sources and sinks. This reflects 
the sensitivity of the correlators and hence the extracted screening masses to the choice of the 
form factor and thus to the particular fuzzing prescription. In the case of maximal smearing 
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L a = 8, L T = 2 the decay of the correlator for tct = (2.9 T) 2 is very close to exponential, 
which is confirmed by the fact that the local masses are nearly constant for 2 < z/a < 6. 
This can be seen from Table 1 and Fig. 5a. However, while the reduction of the error due to 
smearing is impressive for small z, the overestimate is still approximately 12% at z/a = 6. 

Reversing this argument, one could ask for the value of the real screening mass that 
would lead to an observed local mass near the value n{z)a « 1.45. To this end, consider 
first the local masses extracted from the free correlator (Fig. 4c). One observes (cf. Table 
1 and Fig. 5a) that the sensitivity of the local masses to the degree of smearing (and thus 
the fuzzing prescription) is even more pronounced than in the finite mass case, especially for 
small z. However, the striking feature is that, although the theoretically assumed screening 
mass is zero, the local mass extracted via the exponential-type fit ([54]) is finite ! This is 
independent of the degree of smearing, although more strongly smeared sources and sinks 
give smaller values for the local mass. Nevertheless, the local mass extracted via fails 
completely to approximate the real mass if the latter is small compared to T . 

To compare with the results of []IJ we have to bear in mind that our perturbative analysis 
relies on the assumptions that (a) plaquette correlators can be approximated by energy den- 
sity correlators, (b) the main contribution to the energy density correlators is the exchange 
diagram of Fig. Id, and (c) that our smearing prescription accurately models fuzzing as 
employed in [[jj. The limitations of (a) and (b) have already been discussed in Section 2 
and 3. Assumption (c) could be questionable in view of the following fact: the exact shape 
of the source and the sink of the correlator, i.e., the fuzzing or smearing prescription, is 
irrelevant only as long as there is sufficiently large overlap with the wavefunction of the state 
the mass of which is to be measured. This holds in general only for bound states. In our 
case, however, the correlator is supposed to be given by the exchange of two independently 
propagating, dressed gluons. We have already shown (cf. Figs. 4b, c, and 5a) that in this 
case the correlator is very sensitive to the chosen smearing (fuzzing) prescription. 

However, assuming the validity of our assumptions, it is surprising that on the average 
the local masses fi(z) in the case tct = are close to the value obtained in |TJ for the magnetic 
correlator ! In order to get the value of n(z)a ~ 1.45 at z/a = 3 (which was considered to 
give the best estimate for the mass determining the decay of the magnetic correlator in JE]), 
we have to assume ir? to be at most of the order of (2T) 2 (cf. Fig. 5b). On the average, 
n T « T 2 (corresponding to /i rea /a = 0.5) can be best reconciled with the data of j|J. This, 
in turn, means that the local masses may easily overestimate the real masses by at least a 
factor of two. 

However, a more detailed comparison of our results with the lattice data reveals further 
problems with this interpretation. In our calculation, the general behavior of the local masses 
still reflects the remnant enhancement at small z, giving rise to larger fi(z) at small z. This 
trend is not seen in the lattice data and could indicate that our form factor ansatz does not 
accurately model fuzzing as employed in fl| . We suggest that a lattice test of the dependence 
of the local masses on the level of fuzzing could help in clarifying the situation. 

On the other hand, if our fuzzing model works accurately, we are lead to the conclusion 
that either the plaquette correlators may not be approximated by energy density correlators 
or that the latter are not dominated by the two-gluon exchange graph of Fig. Id. If the first 
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possibility is true, then finite size effects are still too large to allow for a simple interpretation 
of the results of |IJ in terms of continuum physics. In the second case, a continuum analysis 
which improves and corrects the simple picture investigated in the present work will be 
extremely difficult, since one does not a priori know which multi-gluon exchange diagrams 
yield the dominant contribution to the energy density correlators. However, in view of 
recent measurements of the spatial string tension jnj] one might speculate that the magnetic 
correlator is dominated by a bound state, as represented by Fig. le. To confirm this, however, 
would require lattice data with much better statistics which is extremely difficult to obtain 
for the plaquette operators studied in U . In any case, if this picture is true, the interpretation 
of the screening mass extracted from the decay of the magnetic correlator as twice the gluon 
magnetic screening mass is not possible. Moreover, this indicates the breakdown of the 



simple perturbative picture of the gluon plasma, in accordance with the conclusions of [£L5 

Finally we comment on the electric correlator. Given the fact that in the continuum 
it is not very different from the magnetic, we may use (]52f) as a rough approximation also 
for the electric correlator. However, in that case the smaller local electric masses found in 
|J (ma ~ 0.8, cf. Fig. 5b), can only be accounted for in our model by assuming that the 
actual electric screening mass is essentially zero! This would be in clear contradiction to the 
screening mass extracted from Polyakov-loop correlations |l| (unless they also suffer from 
ambiguities) . 

We conclude that independent measurements of especially the non-perturbative magnetic 
screening mass are needed to better pin down this physically interesting and phenomenolo- 



gically relevant quantity pD |. 
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Figure Captions: 



Fig. 1: Topologically distinct graphs contributing to the energy density correlator: (a) 
one-loop, two-gluon exchange diagram with free gluon propagators, (b,c) two- and three- 
loop diagrams with free gluon propagators, (d) one-loop, two-gluon exchange diagram with 
dressed gluon propagators, (e) generic ladder diagram representing propagation of a bound 
gluon pair, and (f) one-loop, two-gluon exchange with smeared source and sink. 

Fig. 2: (a) The magnetic correlator in the interaction-free case normalized to T 5 as a 
function of zT (full line) and the contribution of the zero Matsubara mode (dotted line), 
(b) The normalized magnetic (dotted curve) and electric correlator (full line) as functions of 
zT. 

Fig. 3: (a) The magnetic (short-dashed curve) and electric correlator (long-dashed curve) 
for til = (1.4 T) 2 and ttt = (2.9 T) 2 in comparison to the interaction-free result of Fig. 2b. 
(b) Different contributions to the electric correlator: the dotted curve corresponds to the 
contribution of the zero mode, while the full line represents the full result. The exponentially 
damped terms with n ^ constitute the difference between these curves at small zT, the 
anomalous polynomially decaying term the difference for large zT. 

Fig. 4: (a) The magnetic correlator for the interaction-free case in the continuum (long- 
dashed curve: 7r^ = 0, dotted curve: tit = (2.9 T) 2 ) in comparison to a calculation on a 
4 x 16 3 -lattice (squares: m = 0, circles: m = 2.9 T). (b) The magnetic correlator on the 
lattice for m = 2.9 T and different degrees of smearing. The solid curve corresponds to point 
sources and sinks (L a = L T = 0), the dashed curve to smearing with L a = L T = 1, the 
dashed-dotted curve to L a = L T = 2, the long-dashed curve to L a = 4, L T = 2, and the 
dotted curve to L a = 8, L T = 2. (c) The same as in (b) for m = 0. 

Fig. 5: (a) Local masses extracted via (|54j ) from the magnetic correlators of Figs. 4b, c. 
The upper set of curves is for m = 2.9 T, the lower for m — 0. (b) Local masses extracted 
from the magnetic correlator with maximal smearing (L a — 8, L T — 2). The curves are 
labelled with the real gluon mass entering the calculation of the correlator. Black squares 
are local masses extracted from data []J for the magnetic correlator, open circles are local 
masses corresponding to the electric correlator. 
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